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Abstract 

We use the definition of a star (or Moyal or twisted) product to 
give a phasespace definition of the ("-function. This allows us to derive 
new closed expressions for the coefficients of the heat kernel in an 
asymptotic expansion for operators of the form ap 2 + v(q). For the 
particular case of the harmonic oscillator we furthermore find a closed 
form for the Green's function. We also find a relationship between 
star exponentials, path integrals and Wigner functions, which in a 
simple example gives a relation between the star exponential of the 
Chern-Simons action and knot invariants. 



1 Introduction 

The star product on a given symplectic manifold is a very powerful way to 
quantise a given classical theory. Consider a symplectic manifold T, and let 
Aq = C°°(r) be the Poisson-Lie algebra of classical observables. We then 
define a new algebra An which is the quantum analogue of Ao, by letting 

A h = A ®C[[h-\h}} (1) 

and introducing a star product 

f*g = fg + 0{h)J* g -g*f = 0(h) (2) 

and the Moyal bracket 

[f,g}M = f * g - g* f (3) 
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which we demand is a deformation of the Poisson bracket, i.e., 



[f,g] M = ih{f,g}v B + 0{h 2 ) (4) 

The Moyal bracket is also required to be a Lie bracket, thus making A% into 
a Poisson-Lie algebra. It can be proven that twisted products exist on any 
symplectic manifold, and, moreover, that different choices of star product 
corresponds to different choices of operator orderings, |fj]. Furthermore, a 
Woo symmetry exists relating the various choices of star product, ||, and 
each of the various choices specify a characteristic class in the de Rham 
cohomology of V. 

For T = M 2n the unique star product is 

f* g = febH-r^Bg = fg+ l _ lh[f ^ ^ }pB + ^ (5) 

leading to the standard Moyal bracket 

[f, 9\m = 2if sin -}pb^ 9 (6) 

Thus 

,(2«+l)r" A, ,. 4 / n\ 3"/ 3"9 



n=0 k=0 V ' 



k J dq n ~ k dp k dq k dp n - k 

(7) 

The star product comes from the Weyl map which associates a function on 
phasespace to each operator A on L 2 (IR n ), 

A ^ A w {q,p) = J e lup - lvq Tr(U(u,v)A) dudv (8) 

where 

v) = e iup - iv « (9) 

is a translation operator on phasespace giving a ray-representation of the 
Euclidean group E 2n , see 0]. The star product and the Moyal bracket can 
then be obtained from 

A W *B W = (AB) W (10) 
[A W ,B W ) M = (AB — BA) W (11) 
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The inverse of the Weyl map associates an operator / to each function / 
on phasespace 

f w = J e- tup+tvq Il(u,v)f(q, P )dudvdqdp (12) 
and we can then write 

f*9 = (f W 9 W )w (13) 

The Weyl map of the density operator p is referred to as the Wigner function. 
This formalism can be extended to phasespaces which are not just R 2n , see 



2 Star Exponentials and Zeta Functions 

Now, consider a pure state p = we then get the standard formula for 

the Wigner function, W = pw, 

W{q,p) := J e iup - ivq TiU(u,v)pdudv (14) 

{i,{q + y/2)\i,{ q -y/2))e-^^y n (15) 

Let = | A) be a complete set of eigenstates to some Hermitian operator 
H . We can then form the heat kernel of H 

G H (q, q'- a) = (q\e- H °\q') = (q\G H \q') (16) 

and as is well known, Gh determines the zeta function of H, 

Ch(s) = -jy J dacT s - l TrG H (17) 

This provides us with a point of contact with the Wigner- Weyl-Moyal for- 
malism, since, 



TrA = J A w (q,p)dqdp 
(e A ) w (q,p) = Exp(A w (q,p)) 
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where Exp is the *-exponential,0, 



OO OO 

Exp(/) = £ -J™ (18) 



n! s v ' 14 — ' n! 

n=0 „ n=0 



Thus 



( H {s) = -i- / doo* 1 1 iv ' 



r(«) 

= W)S daaS l J E M-Hw{q,p)v)dqdp (19) 

We want to use this formula to find some new expressions for the Schwinger- 
DeWitt asymptotic expansion of the heat kernel. Suppose H = ap 2 + ... we 
can then write 

G H (q, q'; <r) = Y, e ~^ a n(q, q> n ~ d/2 (20) 

n=0 

in d dimensions. This implies 

/OO n n 

G H (q,q;a)dq = ^ / a n (q)a n - d/2 dq = / Exp(-H w (q,p)a)dqdp (21) 
n=0 ^ J 

Now, when = ap 2 + ... then i^vK — a P 2 + ••• where the ellipsis ... stands 
for terms which are at most linear in p. We can then perform the integration 
over p (it is a simple Gaussian integral) to obtain 

J Exp(-(ap 2 + ...)a)dp = (det a)- 1/2 Ti d/2 a- d / 2 E{ q] a) (22) 

where E(q; a) is some function depending on the precise form for H. We 
thus see that the a~ d l 2 factor comes from a Gaussian integration in this 
phasespace formalism. Taylor expanding E(q; a) in a will then give us a q in 
the Schwinger-DeWitt expansion. Notice that such asymptotic expansions 
are rather easily derived with this formalism. For instance, if H = ap 1 + v(q) 
then we get a factor a~ d l 1 for 7 > 1 from the integration over p, and we 
can then find a generalised Schwinger-DeWitt expansion for such (pseudo- 
differential operators. 
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Also interesting is an expansion in powers of h. Suppose / is independent of 
h, we then write 

00 

Exp(/) = J2 h2n ^n(f) (23) 

n=0 

By writing down the definition of the twisted product, we see that only even 
powers of H can occur - odd powers only become relevant if the exponent 
itself is ^-dependent. Furthermore, the h 2 term, say, will receive contributions 
form all the twisted powers of /. More precisely, we see 

r = f n__J2 /^ 2 (/, /""I"*) + (h A ) (24) 
fc=0 

The higher orders of H will not give so simple contributions. The fourth 
power, for example, will receive contributions not only from f k oo4(f, f n ^ 1 ^ k ^ 
analogous to the /i 2 -terms above, but also from uj 2 (f k , oj 2 (f\ f n ~ k ~ 1 )) and 

, , / fk fl\, , I fm fn—m—k—l\ 
U2\J , J )^2\J , J )■ 

From this it follows that we can write 

00 -, n— 2 



w) = -^E^E/W/,/"- 1 -*) 



8 ^ n\ 

n=2 k=0 

= -\^(f, f)F 2 (f) - \u2{f, f)G 2 (f) (25) 

where we have used 

u 2 (f, f m ) = mr-'^if, f) + m(m - l)f m ~ 2 u 2 (f, f) (26) 

with 

~ (f f) ._ d2 f ( d f\ 2 2 d2 f d f d f 92 f ( d f\ (27) 
dq 2 \dp J dqdp dq dp dp 2 \dq J 

The functions F 2 , G 2 turn out to be 

00 n— 2 00 

w) - z^("-i-*)/- 2 45:^/"- 2 <»> 

n=2 k=0 n=2 v ' 
°° 1 n ~ 2 1 °° Q 

G >(/> = £^5> .-!-*)(»- 2 -*)/-» = i^.»z» /-tfc) 

n=2 ' k=0 n=2 V >' 
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One will be able to write in a similar fashion 

w) = ^ m/, m(m(f, f)G,(f) + uj 2 (f, f) 2 F 2A (f)+ 

f) 2 G 2A (f) + /)*M/)) (30) 

for suitable functions F 4 , G 4 , F 2i4 , G 2i4 , i?2,4- Clearly, as the power of h grows, 
these expressions become more and more cumbersome. If one would wish to 
go to arbitrary high powers of H it would be worthwhile to find a diagram- 
matic expression and some "Feynman rules" for constructing £ n . It is clear 
from the recursive way these are constructed that such "Feynman rules" must 
exist. We will not attempt to find these rules here, however, but merely re- 
strict ourselves to noting that often one can assume that powers of / vanish 
for sufficiently high powers, or at least become very simple. An example of 
this could be / = ap 2 + v(q), if the potential is not too big, one can write 
jn ^ a np2n f or n sufficiently large. Since this particular example is of great 
physical importance - it gives the Hamiltonian of a Schrodinger particle for 
a = I /2m or the Klein-Gordon field coupled to some external field in a non- 
minimal way (e.g. curvature) with a = 1/2, v{q) = -|- + £R(q) + V(q) - we 
will write down the explicit result for £ 2 . From the general formula we see 

£ 2 (ap 2 + v(q)) = -\av"F 2 - -(2a 2 v"p 2 + av' 2 )G 2 (31) 

The next term, £4, will contain third and fourth derivative of v(q), hence for 
v ~ q k for some positive integer k the result will simplify since the derivatives 
of v will vanish sooner or later. For k = 2, the harmonic oscillator, v(q) = 
\mujq 2 , a = l/2m we have already 

£ 2 = --uF 2 - \ (-p 2 + mu 2 q 2 ) G 2 (32) 
4 8 Vm / 

We have plotted e* + £ 2 {f) for / = p 2 + q 2 , h = 1 in figure 1. 
Another interesting case is v(q) — z/q corresponding to the Coulomb poten- 
tial (again a = l/2m). Here £ 2 becomes 

and + £ 2 (f) has been plotted for / = p 2 + 1/q in figure 2. 
A Yukawa-like potential v(q) = z exp(— yug)/g would lead to 

£2 = -^e-^ + ^ + ^V 1 )^- 
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^-^e"^ [2p 2 (q + m 2 + - A V + mz(l + 2 m + ,V) ) G 2 

We have not plotted this, since the exponential factor will only dampen the 
behaviour a bit but essentially it would look the same as for fi = 0. 
One should note that for very large values of n, the sums defining F 2 , G 2 will 
simplify to essentially J2 n f n / n ^ thus 

F 2 (-f) = e-^ 2 (-f) G 2 (-f) = e-t&i-f) (34) 

where J r 2 , Q 2 have slower growth than the exponential, making F 2 , G 2 inte- 
grable as functions of p. This will be needed for the relationship with the 
Schwinger-DeWitt expansion to make sense. Moreover, for / = ap 2 + v(q) 
the integral over p will give a~ d l 2 times a function of q by simple Gaussian 
integration. Once more we see that the factor o~ d l 2 which "shifts" the pow- 
ers of a in an asymptotic expansion of the heat kernel away from the simple 
Taylor series arises in a straightforward manner from the phasespace formal- 
ism upon integrating out the momentum variable. 

A way to systematically compute such "quantum corrections" can be found 
by noting that the star exponential is a solution to the differential equation 

§ = (35) 

the solution of which is precisely G = Exp(-fa). Now similarly g = e~f a 
is the solution of a differential equation obtained by replacing the twisted 
product in the above equation by an ordinary one. Write 

G(q, p; a) = Exp(-/(g, p)a) = g(q, p; a)G(q, p; a) = e~^>g{q, p; a) (36) 

then we have that Q is a solution to 

^ = g- 1 f*{gQ)-fQ (37) 

= -^hg- 1 {f,gg} PB + ^g- 1 LU 2 (f,gQ) + ... (38) 

where we have used the definition of the twisted product. Since g = e~^ a the 
factor g can be pulled out of the Poisson bracket cancelling the g~ l factor 
outside. For / = ap 2 + v (q) we get 

2 '2 o/„ „„2 2\ 



ou 2 (f, g g) = 2ag \{a 2 v 2 - 2(a - ap 2 u 2 )v")Q- 

2pav"^- + — v"^- - 2av'^- + ^ 
dp 2a dp 2 dq dq 2 



(39) 



7 



Supposing Q can be Taylor expanded in powers of h, 

CO 

g = Y^g n h n (40) 

n=0 

with Q Q = 1 we get the following recursive relation 



n 



-^Gn = ~J2 ^Fk} 9 W/, SGn-k) (41) 
k=l 

from which we see that Q\ = as expected and that Q 2 becomes simply 

Q2 = ~\ f ef°u 2 {f,e-f°)da (42) 
which for the chosen / can be integrated quite readily to give 

g 2 = QcktV - ^aa 3 (v' 2 + 2apV)) (43) 
By construction we also have 

02 = e fa S 2 (44) 

leading finally to 

£ 2 = \-e- v °^a z l 2 {2v" - av' 2 ) (45) 

In a similar way, one can find analytical expressions of the remaining £ n 's. 
In any case, we can use the expression for the S 2n to find expressions for the 
Schwinger-DeWitt coefficients order by order in h. Examples of this will be 
given later. 

For completeness we will quickly list the formula for Q±. Straightforward 
computations yield 

Qa = ^^(/.e^^J + ^eNae^))^ (46) 

= l a V 3 t,W(5 - 15apV + Aa 2 p A a 2 ) - —a 2 a\v 2 + 2a P 2 v") 2 + 
480 V y F J 288 V y ' 



1 

-a 



V(9vV + 18ap 2 v" + Aap 2 v'v^) + -UW 4 ) - 



240 ' 48 

— a 2 a 2 (5v" 2 + Av'v {3) + ap 2 v (4) ) (47) 
96 



From this one can then compute Ei n as 



£ 2 n = e~ fa Q2n (48) 

As will be clear from these simple examples, the computations are all rather 
elementary, allowing one rather quickly to find all the relevant f^n's, and, we 
will see, also the coefficients in a Schwinger-DeWitt asymptotic expansion 
for the heat kernel. 



3 Path Integrals, Wigner Functions and Lo- 
calisation 

We can make one more important connection. Let \() be a set of coherent 
states (i.e., for the Heisenberg algebra in n dimensions, h n , we have ( G C n ). 
We then have 

Exp(-H w a) = J fiCC.Oe-^^^C.O ( 4 9) 

where 

0(C,C'):=^^ (50) 

s is a parameter along the path in £-space and T>((, (') is the functional 
measure 

v(c, c)--= iim n^w^r 1 ^ 

with Co = (, (n = ('■ See for instance 0. This formula can be read two ways, 
either the functional integrals allows one to compute the star exponential, or 
the star exponential allows one to define a regularised functional integral. 
Introducing the eigenfunctions |A) of H, we can write H = 52 A A]A)(A|, but 
we can also write the heat kernel as 

GH = e-^ = 5>- ACT |A)(A| (52) 



Thus we arrive at 



{e- H °) w = Ex V (-H w (q,p)a) (53) 



A 
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where we have defined 

W X =(\X)(X\) W (55) 

This allows us to "localise" a functional integral, turning it into a sum over 
a discrete set (the spectrum of the operator). It also says that the Weyl 
transform of the heat kernel is a kind of Laplace transform (A — > er) of the 
Wigner function. 

From the heat kernel one can also compute the Green's function G = H^ 1 
by simply integrating over a, 

POO 

G = H- 1 = - / e' Ha da (56) 
Jo 

and this implies that we can express the Green's function in terms of the 
Wigner function as (using the linearity of the Weyl transform) 

G(q,p) = - J>-Hy A ( ? ,p) = (H-^w := H^ 1 (57) 

A 

or 

/■oo 

G(q,p) = Ex P (-H w (q,p)a)da (58) 
Jo 

giving us an interpretation of the integral of the star exponential. 



4 Applications 

Our first application of the relationship between path-integrals and star- 
exponentials involves quantum mechanics. Let H = ap 2 + v(q) be the Hamil- 
tonian for a Schrodinger particle in d = 1 dimension. Then 

U(u,v) = e iu P- iV(i (59) 

is the standard translation operator, giving a ray representation of the group 
of translations in two dimensions (i.e. on phasespace), E 2 . We then have 

Exp(-(ap 2 + v{q))a) = j e^-^' e~ ^ a ^ +v ^' )ds V(C,X') (60) 

= e -(°P 2 +«(*)V + h 2 S 2 (Ha)+0(h 4 ) (61) 

A 
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here one can interpret a as an inverse temperature, a = (3. In this case, 
the eigenvalues A are the energies and the W\ become the Wigner functions 
of the eigenstates, i.e., the solutions to the time-independent Schrodinger 
equation. For the explicit example of an harmonic oscillator we then arrive 
at the relations 

E X p(-I(^ + WgV) = y"e^-^'e-/o CT ^+TC' 2 ^p(C,C') (63) 

/r / ^ + m v)G 2 + 0(fi 4 ) (64) 
m J 

Y,W n (p,q)e-^ n+1 ^ (65) 



n=0 



W n M = ^re -^--V^f^f^ + ^V)) (66) 



where 

fmn\ n \ fvjj \m 

is the Wigner function for the harmonic oscillator, ||. This also provides 
us with an explicit formula for the star exponential and hence of the path 
integral. The expression in terms of the Wigner function (i.e., in terms of 
Laguerre polynomials L n ) show that we can "localise" the functional integral 
on a countable set, in the same way one can "localise" a curve integral in the 
complex plane on a finite set of poles, the value of the integral being propor- 
tional to the sum of residues at these poles. Here we have an integral over 
a set of more than continuum cardinality (the set of paths in the complex 
plane) instead of just a continuum (a specific curve in the plane), hence the 
set on which we "localise" is no longer finite but countable. In general one 
would have the functional integral localised on a set of at most continuum 
cardinality (corresponding to a spectrum with continuous eigenvalues and 
not just discrete one). 

The above relation between the star exponential and the Laguerre polyno- 
mials also imply that the Green's function for the harmonic oscillator can be 
written as a mode sum 



G{q,p ) = 2-^—pr- ( — , -, /0 n , e V sL n [ — [— + mwq 

irn z uj{n + l/2)n\ \nui \m 



(67) 



11 



This is then the phasespace Green's function of this important Hamiltonian. 
It turns out that we can find a closed formula for this sum. The generating 
function for the Laguerre polynomials is 



71=0 



Consequently, 

00 ( — -\\n r +1/2 



n=0 

from which we can see at once that 



G(q,p) = —£(x)e-* x (69) 



\e- x ' 2 lim f" 3 ' 2 [ ^e-^dt (70) 

nh *— 1 y 1 - * v ; 



with x = j|j + mw 2 g 2 j . Unfortunately, I have not been able to compute 

this analytically but it certainly straightforward to do so numerically. 
We can also find the Schwinger-DeWitt coefficients a n in this example. By 
performing the integral over p and Taylor expanding the result in a we get 
at once 

a = V2^ + 0(h 2 ) (71) 
Gl = -nm 2 ujq 2 + 0(k 2 ) (72) 

a 2 = ^vrmVV + 0(fr 2 ) (73) 

and so on. 

For an arbitrary potential in one dimension we would get 

a = V2^+0(h 2 ) (74) 
ai = -y/2irmv(q) + 0(h 2 ) (75) 
a 2 = ^V + 0(/l 2 ) (76) 

The h 2 terms will contain derivatives of v and we see that to lowest order in 
h we can treat the potential as a constant. In all cases, the 0(J% 2 ) terms are 
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given by f S 2 dp. 

Inserting our expression for S 2 found above, for ap 2 + v(q), we can readily 
compute the momentum integral 

J E 2 dp = ^V^(r 3/2 e~ av (2v" -av' 2 ) (77) 

Notice that the following powers of o will appear: 3/2,5/2,7/2, Thus 

the power a^ 1 ^ 2 appearing in both the phasespace Taylor series and in the 
asymptotic Schwinger-DeWitt expansion comes, as was to be expected, solely 
from the "leading symbol" p 2 , the remaining terms only giving higher powers 
of a. 

The 0(h 2 ) correction, S 2 a n , to the Schwinger-DeWitt coefficients then be- 
come 

5 2 a = 5 2 ai = 5 2 a 2 = -y/anv" (78) 

Similarly, by integrating £ 4 one can find the 0(ft 4 ) corrections, which will 
contain even higher powers of a. 

One of the problems with the standard asymptotic expansion is that it ig- 
nores boundary conditions and other global properties of the spacetime man- 
ifold and of the fields. This is a problem, for instance, in computations in 
curved spacetimes where one needs some such global contribution (coming 
from the vacuum definition) when one wants to find the renormalised energy- 
momentum tensor. While definitely not proven, there is some hope that the 
present formalism can cure some of these problems, this is so because the 
expressions for the heat kernel involves an integral over momentum space, 
i.e., on a curved manifold the cotangent bundle, and this has some infor- 
mation stored in it about the global properties of the spacetime manifold 
(how the bundle is glued together for instance). A full investigation of this 
problem lies beyond the scope of the present article which only claims to lay 
the foundation and to point out the applicability of the phasespace formalism. 

As another simple example we will consider where the canonical 

variables are ir l a ,A?, with % = 1,2,3 a spatial index and a = l,...,n a Lie 
algebra index, n = dimjj. Formally, then, the IT-function is given by 

Yl — e i{u,A)-i(v,ir) (^jgj 

where u G (VL l <g> g)*,t> G VL l (g) g are objects "dual" to A,n. The brackets 
denote the pairing between the spaces of {A, n} and their duals. We can 
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thus write 

(u,A) = J uiAfdx (v,n) = J vtKdx 

We will be interested in the particular case where 

u\dx = d^\ v^dx = dr)i (80) 

where then 7, t] are loops. In this case II becomes the phasespace analogue of 
a Wilson loop. As the observable H we will pick the Chern-Simons functional 

H = -ig J tx{A A dA + A A A A)dx = -iScs[A] (81) 

We then have 

(U) t := J e if i A - i *>i v e- 4 fi Sos[A]dt VAThr (82) 

This path integral describes a family (indexed by t) of Chern-Simons theories. 
Writing 

^ n = J Tr7r ■ f]ds = J Tin ■ Adx (83) 

where A° = 5(x,r](s))rjf is the so-called form factor of the loop, we can 
perform the ir integral obtaining 

(U) t = S(A(rt)) J e l ^ A - iJ o Scs[A]dt VA (84) 

Two particular cases are interesting, one is where {A t } is a constant family, 
i.e., A t = A = A, Vt, the other is where it is a family centred on a finite 
number of points, i.e., A t 7^ only for t = t , ...,t N . In the first case, we have 
Jq Scsdt = tScs so putting k = Angt and q = exp(|||) we have, according 
to Witten, [|, 

ExpHS cs ) = (n) = S(A( V ))c(k)- w ^ J q {j) (85) 

where 1^(7) is the writhe of the loop, and J q is a Jones polynomial. This 
also imposes the quantisation condition k = Angt G Z. The second case will 
give a finite sum over such expressions with different coupling constants kt 
and consequently indices qt. We can consider the general expression as an 
interpolation between knot invariants for different values of k, q. 
The equation (|85|) gives us an interpretation of the star exponential in this 
instance, and potentially also another way of computing Jones polynomials 
and perhaps of finding relations between them. 
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5 Conclusion 



We found a relationship between zeta functions, path integrals, star expo- 
nentials and Wigner functions. This allowed us to regularise path integrals, 
derive phasespace expressions for determinants and effective actions and fi- 
nally to localise path integrals in the sense of topological field theory (i.e., 
write a functional integral as a sum over a set of lower cardinality - either 
discrete or continuous). By Taylor expanding the phasespace expression for 
the heat kernel, we arrived at a comparison with the standard Schwinger- 
DeWitt expansion. We computed the Schwinger-DeWitt coefficients of the 
heat kernel to 0(h 4 ) explicitly for / = ap 2 + v(q). Another explicit example 
was taken from topological field theory, where we constructed the star expo- 
nential of the Chern-Simons functional in three dimensions, thereby getting 
the Jones polynomials. The application to the heat kernel, however, seems 
to be the most promising. 
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Figure 1: The star exponential, Exp(/) to 0(ft 4 ), for / = p 2 + q 2 in units 
with h = 1. 



Figure 2: The star exponential for / = p 2 + l/q to 0(h 4 ). 
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